Ab Initio Study of ^°Ca with an Importance Truncated No-Core Shell Model 
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We propose an importance truncation sciieme for the no-core shell model, which enables converged calcu- 
lations for nuclei well beyond the p-shell. It is based on an a priori measure for the importance of individual 
basis states constructed by means of many-body perturbation theory. Only the physically relevant states of the 
no-core model space are considered, which leads to a dramatic reduction of the basis dimension. We analyze 
the validity and efficiency of this truncation scheme using different realistic nucleon-nucleon interactions and 
compare to conventional no-core shell model calculations for ''He and '''O. Then, we present the first con- 
verged calculations for the ground state of '"^Ca within no-core model spaces including up to 16ftfi-excitations 
using realistic low-momentum interactions. The scheme is universal and can be easily applied to other quantum 
many-body problems. 

PACS numbers: 21.60.Cs, 21.30.Fe, 21.60.-n 
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The theoretical description of nuclei as strongly interact- 
ing quantum many -body systems bears two fundamental chal- 
lenges: (i) The interaction between the constituents is a com- 
plicated, residual force resulting from the strong interaction 
of quantum chromo dynamics (QCD). (ii) The nuclear many- 
body problem, which is dominated by correlations, has to be 
solved reliably for systems containing from 2 to a few hun- 
dred nucleons. There has been substantial progress on both 
aspects over recent years. In addition to traditional realistic 
nucleon-nucleon interactions, like the Argonne V18 ill] and 
the CD Bonn potential |01, which reproduce the experimen- 
tal nucleon-nucleon scattering data with high precision, in- 
teractions derived in the framework of a chiral effective field 
theory became available |[3),|4J. The latter provide a system- 
atic and consistent picture of two- and three-nucleon interac- 
tions based on the symmetries of QCD. Given these interac- 
tions the solution of the nuclear many-body problem remains 
a formidable task. For light nuclei, ab initio methods like 
the no-core shell model ISila] or the Green's Function Monte 
Carlo method 10, Ht] provide tools to solve the many-body 
problem numerically. Recently, the former has been used to 
study the properties of mid-p-shell nuclei based on chiral two- 
plus three-nucleon interactions for the first time |@]. 

Methods like the full no-core shell model (NCSM) become 
computationally intractable for heavier nuclei. Currently con- 
verged NCSM calculations are limited to the p-shell. For 
heavier nuclei one relies on severe approximations, such as 
perturbative expansions on top of a Hartree-Fock solution 
i Tol n], or completely phenomenological approaches. Only 
few approaches, e.g. the coupled-cluster method Iil2,il3ll. are 
able connect the domain of strict ab initio methods with the 
domain of approximate models for heavier systems. 

In order to extend the scope of the ab initio NCSM to heav- 
ier nuclei we propose an extension using a priori information 
on the importance of individual basis states obtained from 
many-body perturbation theory. We use this importance mea- 



sure to truncate the NCSM model space to the physically rel- 
evant basis states and thus reduce its dimension substantially. 
The current limitation of full NCSM calculations results from 
the sheer dimension of the matrix eigenvalue problem. Typ- 
ically the many-body basis, which is represented by Slater 
determinants of harmonic oscillator single-particle states, is 
truncated at a maximum number iV„iax of harmonic-oscillator 
excitation quanta thus defining the so-called iVmaxfif^ model 
space. Its (m-scheme) dimension grows exponentially with 
the number of nucleons A and the truncation level A^max- 
For ^^O this limits the presently tractable model space to 
^max = 8 lfl4ll corresponding to an effective m-scheme di- 
mension of 6 X 10®. For *"Ca the dimension of the S/iil 
model space is 2 x 10^^ — well beyond the capabilities of cur- 
rent shell model codes. Many of these basis states are irrele- 
vant for the description of any particular eigenstate, e.g. the 
ground state. Therefore, if one were able to identify the im- 
portant basis states beforehand, one could reduce the dimen- 
sion of the matrix eigenvalue problem without loosing predic- 
tive power. This, as we propose in this Letter, can be done 
using many-body perturbation theory. This type of problem 
is not restricted to nuclear physics but common to many ar- 
eas of quantum many-body theory. For the description of 
strongly corelated lattice systems, e.g. in the context of ul- 



tracold atomic gases in optical lattices 111511 . or other atomic 
or condensed matter systems, one encounters the same limita- 
tions of exact diagonalization approaches. 

We start out with a translationally invariant Hamiltonian 
Hint composed of the intrinsic kinetic energy Tint = T — Tcm 
and a two-nucleon interaction Vnn- In the first part of this 
Letter, we use the phase-shift equivalent correlated interac- 
tion VucoM discussed in Refs. ifioi For three- and 
four-nucleon systems it provides binding energies in good 
agreement with experiment without the inclusion of an ex- 
plicit three-nucleon force. Together with its good conver- 
gence properties it provides a suitable benchmark interaction 
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for the present study. In the second part we use the Mowfc low- 
momentum interaction lll7ll to discuss the ground state struc- 
ture of and "^^Ca. 

The general concept of the importance truncation is as fol- 
lows: We start with a reference state l^'i-ef), which represents 
a zeroth-order approximation of the many-body state we are 
interested in. Here, we limit ourselves to the ground state, but 
the scheme can easily be extended to excited states or more 
than one reference state. Starting from l^'ief) we can build a 
many-body space by generating all possible n-particle-7i-hole 
(npnh) excitations up to the excitation energy A^,„ax^i^- By 
increasing n {< A), we eventually recover the translationally 
invariant Nan^^hn model-space of the NCSM. We now esti- 
mate the contribution of a given npnh-state to the ex- 
act eigenstate after the NCSM diagonalization via many-body 
perturbation theory. In first-order MBPT the amplitude of the 
state 1$^) is given by 



0.02 



(1) 



where H' is the Hamiltonian of the perturbation and e^^ are 
the unperturbed energies of the two configurations. When 
using the harmonic oscillator basis, the unperturbed Hamil- 
tonian is just the one-body harmonic oscillator Hamiltonian, 
Hq = Hho, and the perturbation is H' = Hint — Hho, where 
the contribution of Hho to the matrix element in ([TJ vanishes. 

If we restrict ourselves to two-body interactions, then H' 
contains only one- and two-body terms such that van- 
ishes for 3p3h and higher-order configurations. In principle 
higher orders of perturbation theory are required to gener- 
ate states beyond the 2p2h level directly. Since this becomes 
computationally inefficient, we resort to an iterative scheme. 
In a first iteration we generate Iplh and 2p2h states start- 
ing from the reference state, retain those with an importance 
weight Ki, > Kmin, and solve the eigenvalue problem in this 
space. In the next iteration, the dominant components of the 
ground state |5'o) = 51! ^ C*!/ \^l>) obtained from the diago- 
nahzation in the previous step are used as reference state, i.e. 

I*ref) = EI'^"'"^'" 1*!^) with a-ef « 0.005. Since this 
state already contains up to 2p2h admixtures one obtains non- 
vanishing importance weights ([T]i for states up to the 4p4h 
level. After applying the importance truncation > Knun 
the eigenvalue problem is solved in the extended space. This 
cycle can be repeated until the full Nnm^fiil model-space is 
generated in the limit Kniin = . For the first applications 
presented here, we restrict ourselves to the 4p4h level, exten- 
sions will be discussed elsewhere. 

The correlation plots presented in Fig. [T] demonstrate that 
this prescription provides a reliable a priori estimate for the 
size of the a posteriori amplitudes of the individual ba- 
sis states after solving the eigenvalue problem. There is a 
clear correlation between the coefficients and the impor- 
tance measure k^. Therefore, the additional truncation of 
the many-body space to configurations with k^, > Kniin pro- 
vides an efficient means to remove irrelevant configurations 
from the outset. Eventually the eigenvalue problem within 
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FIG. 1: (color online) Correlation between the a priori importance 
weight and the a posteriori amplitudes Cu resulting from the ex- 
plicit diagonalization. Shown are the 2p2h (a) and 4p4h contributions 
(b) for a 6hQ, calculation for ^"^O with Vucom at hO, = 20MeV. The 
dashed vertical line marks a typical truncation parameter Umin- 
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FIG. 2: (color online) Dependence of the ground-state energy of 
'^^O (a) and '"'Ca (b) on the truncation parameter Kmin obtained with 
VucoM. The three data sets correspond to model spaces with up to 
2p2h (•), 3p3h (♦), and 4p4h states (■), respectively. The lines show 
a 5th order polynomial interpolation. 



the truncated space is solved without further approximation. 
Hence, formal properties such as the variational principle and 
the Hylleraas-Undheim theorem remain intact — in contrast to 



other approaches like the coupled cluster method 11121 11311 . 
When reducing the truncation parameter Kj^m, we are guar- 
anteed to converge from above to the exact eigenvalues for 
the given interaction. 

The dependence of the ground state energy on the trunca- 
tion parameter Kmin is illustrated in Fig. |2]for ^^O and ^°Ca. 
The different data sets correspond to calculations with spaces 
comprising up to 2p2h, 3p3h, and 4p4h states, respectively. 
We always observe a very regular behavior, which allows for 
an extrapolation of the energies towards Kmin — correspond- 
ing to the complete model space. We use a 5th order polyno- 
mial fit that turns out to be very robust. 

A first indication of the efficiency of the importance trun- 
cation scheme is given in Fig. |3] Shown are the ground state 
energies of "'He and ^^O as function of iVmax computed for 
the Vucom interaction using the importance truncation (after 
extrapolation to Kmin = 0) in comparison to the full NCSM 
calculation performed with the Antoine code llSll . For "'He 
the truncated calculation including up to 4p4h configurations 
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FIG. 3: (color online) Convergence of the ground-state energy for 
*He (a) and ^"^O (b) versus A'max obtained using Vucom- Shown 
are three data sets corresponding to model spaces with up to 2p2h 
(•), 3p3h (♦), and 4p4h-states (■), respectively. Black crosses (+) 
indicate the results of full NCSM calculations. Lines to guide the 
eye. 

generates the full iV,„ax^i^ model space in the limit K^in — 0. 
The excellent agreement between full and truncated calcula- 
tion in this case proves the reliability of the importance trun- 
cation scheme. Nonetheless, the dimension of the eigenvalue 
problem is drastically reduced by the importance truncation: 
For A^max = 20 the dimension of the full NCSM model space 
is beyond 2 x 10' whereas the dimension of the truncated ba- 
sis is of the order 10^. This redcution becomes even more 
striking for heavier systems. 

For ^^O a similar picture emerges, the importance truncated 
calculations up to the 4p4h level are in very good agreement 
with the full NCSM. For the 8hi} model space the full calcula- 
tion yields a slightly lower ground-state energy, which is due 
to states beyond the 4p4h level. The relevance of 4p4h corre- 
lations is a specific property of ^^O (a-clustering) and reveals 
through the compareably strong effect of those states on the 
ground-state energy. 

The A'^max-dependence of the ground-state energy of ''"Ca 
for two different oscillator frequencies is presented in Fig. ID 
In comparison to ^^O, the impact of 4p4h configurations is 
negligible and it is expected that higher-order configurations 
will not affect the ground-state energy. For the oscillator fre- 
quency hfl = 17 MeV even the inclusion of 3p3h states does 
not lead to a significant effect. An exponential extrapolation 
of the 3p3h energies for hft = 17 MeV to infinite model- 
space size gives Eoo ~ —316 MeV, which can be compared 
to the experimental binding energy of —342.05 MeV. Keep- 
ing in mind that our calculation — ^because of the restriction to 
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FIG. 4: (color online) Convergence of the ground-state energy of 
*°Ca as function of model-space size A^max for HQ = 17 MeV 
(lower curves, left-hand axis) and 20 MeV (upper curves, right-hand 
axis) using the Vucom interaction. Symbols as described in Fig. [3] 

the 3p3h level — provides an upper bound, this result is very 
encouraging. It shows that the phase-shift equivalent two- 
nucleon interaction Vucom is able to provide a realistic de- 
scription of binding energies also for heavier nuclei without 
the inclusion of an additional three-body interaction. 

In this intermediate mass regime the importance truncation 
works very efficiently: The dimension of the importance trun- 
cated space used for ^"Ca around iV,nax = 16 is of the or- 
der of 10^, which is smaller than the full NCSM space at 
^max = 4. In order to warrant that contributions of spu- 
rious center-of-mass excitations to the energy are negligible, 
we monitor the expectation value of the operator H^™ — |ftfi. 
We obtain typical values of 100 keV and below proving that 
the importance truncation does not generate spurious center- 
of-mass excitations, which are otherwise absent in a complete 
NniiLxfi^ model space. 

Following these benchmark calculations with Vucom, we 
apply the importance truncation for a systematic study of the 
ground states of ^^O and '"'Ca with the Vfowfe low-momentum 
interaction derived from the Argonne V18 potential with a 
cutoff momentum A = 2.1 fm~^ lfl9ll . The T^owfc interaction 
shows an even faster convergence than Vucom and therefore 
provides a useful reference for a quantitative comparison of 
different many-body approaches. 

The convergence of the ground-state energy of ^^O and 
^''Ca as function of the oscillator frequency hfl for different 
model-space sizes from iVmax = 4 to 14 is illustrated in Fig. 
[5] por i^o configurations up to the 4p4h level have been in- 
cluded, for ^°Ca up to 3p3h configurations are considered. 
For both nuclei we observe a very rapid convergence of the 
energies, such that fully converged results are obtained. Ta- 
ble H] summarizes several ground state observables computed 
for A^max = 14 and oscillator frequencies around the energy 
minimum. In addition to the ground-state energy E for this 
model space, the energy Eoo resulting from an exponential 
extrapolation iVmax ^ oo is given. For ^^O this result can 
be compared with a recent coupled cluster calculation using 
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FIG. 5: (color online) Ground-state energy of 0(a)and Ca(b)as 
function of the harmonic oscillator frequency HQ for different model- 
space sizes up to A^max ~ 16 obtained with the Viowfc interaction. 
For up to 4p4h configurations are included, for ''"Ca up to 3p3h. 
Crosses indicate the results of full NCSM calculations. 

the same interaction, which yields an extrapolated energy of 
-142.8 MeV III. This is in good agreement with our result 
keeping in mind that our calculation provides a variational up- 
per bound, since configurations beyond 4p4h are not included. 
Furthermore, Tab. H] lists the translationally invariant root- 
mean-square radii (point nucleons) as well as the charge radii 
(including proton and neutron form factors). From the eigen- 
vectors in the importance truncated space we have also ex- 
tracted charge form factors and translationally invariant den- 



TABLE I: Ground-state energy E, point mass radius iiims and charge 
radius obtained for ^"^O (up to 4p4h for iV^ax = 14) and '"'Ca 
(up to 3p3h for A'^max = 16) using the Viowfc interaction for dif- 
ferent oscillator frequencies. Eoo indicates the ground-state energy 
obtained from an exponential extrapolation. 
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sity profiles, which provide an independent means for extract- 
ing translationally invariant radii. In comparison with exper- 
iment, the two-body interaction VJowfc overestimates the bind- 
ing energies and underestimates the radii significantly. 

In conclusion, we have presented an importance trunca- 
tion scheme which reduces the model space of the NCSM 
to the physically relevant states. The relevance of individual 
basis states is quantified using an a priori importance mea- 
sure derived from many-body perturbation theory. This novel 
scheme has proven very robust and extends the range of re- 
liable NCSM calculations to regions of the nuclear chart far 
beyond the p-shell. In this Letter, we have shown the first con- 
verged NCSM calculations for the ground state of *"Ca with 
two different realistic NN interactions. The scheme is univer- 
sal and can easily be adapted to other quantum many-body 
problems treated by diagonalization techniques. 
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